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We show by microscopic calculation that thermodynamics of the multicomponent Sutherland 
model is equivalent to that of a free particle system with fractional exclusion statistics at all temper- 
atures. The parameters for exclusion statistics are given by the strength of the repulsive interaction, 
and have both intra- and inter-species components. We also show that low temperature properties 
of the system are described in terms of free fractional particles without the statistical parameters 
for different species. The effective exclusion statistics for intra-species at low temperatures depend 
on polarization of the system. 



Q\ . I. INTRODUCTION 



State counting rule of a large number of identical particles depends on the quantum statistics; A number of bosons 
can occupy a single particle state ( no exclusion ) while the presence of a fermion in a state excludes other fermions 
(complete exclusion) . Although particles and conventional elementary excitations obey bosonic or fermionic statistics 
so far, elementary excitations with fractional statistics attracted theoretical interests recently. Stimulated by the 
fractional quantum Hall systems and one-dimensional quantum spin chain, Haldane proposed a new type of quantum 
statistics from the viewpoint of state-counting. Q He defined a fractional generalization of the Pauli's exclusion 
, principle so that it gives a state-counting rule interpolating the bosonic and fermionic one. This new type of quantum 
statistics is called the fractional exclusion statistics (FES). Subsequently, thermodynamics of single component free 
particles which obey FES was constructed, based on the new state-counting rule. Thermodynamics was derived 
for the system of multicomponent particles obeying FES in ref. 4. All the above works start from the state-counting 
rule or entropy. Hence in refs. 1-4, the question how the FES is realized in concrete models is left open. 

The validity of FES was tested in two solvable models defined in the continuum space: one-dimensional bose gas 
interacting with the delta function and the Calogero-Sutherland model In those papers, it was found that the 

thermodynamics of the two systems are described in terms of free particles with FES. The above models, however, 
have no internal structures and hence the following question arises: How the FES is realized in concrete models with 
internal structures? 

In this paper, we establish by microscopic calculation the relation between FES and multicomponent Sutherland 
model (MSM). We show that the thermodynamics of MSM is the same as that of free particles with FES at all 
temperatures. Furthermore we obtain a simpler description with FES in low temperature region. 

This paper is organized as follows. In the next section, we review FES and introduce the statistical parameters. In 
section 3, we construct the thermodynamics of MSM by following the Yang- Yang's method. In section 4, we show that 
the MSM is thermodynamically equivalent to a free particle system with the mutual statistics at all temperatures. 
Furthermore, in section 5 we find that the MSM reduces to a free particle system with fractional statistics but without 
jjj ■ the mutual statistics (which is called "g-on " ||). 
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II. FRACTIONAL EXCLUSION STATISTICS AND THERMODYNAMICS 



General formulation and thermodynamic properties of FES were presented by several authors. |l]|3| In this section, 
we give a brief review on the FES, in order to fix the notation and definition of terminology. First divide the total 
system into many macroscopic subsystems each of which is labeled by a. Here we denote by D a the number of available 
single particle states in the subsystem a in the presence of other particles, and by N a the number of particles in the 
subsystem a. Haldane (!]] introduced the statistical parameter g a p through the following differential relation: 

AD a = -J29*f3&N . (2.1) 



Integrating eq. (2.1) with N a , we obtain 



D a = G a -y2g aP Np, (2.2) 
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where G a is the number of single particle state in the absence of other particles. The total number of microscopic 
states with {N a } is given by 



11 NJ.DJ. ■ 



(2.3) 



If we set 9a@_= 0, equation ( |2.3| ) reduces to the boson statistics in the thermodynamic limit. In t he c ase of g a p = 8 a p, 
equation ( |2.3| ) reduces to the fermion one in the thermodynamic limit. Here we can regard eq. (2.3) as the fractional 
generalization of conventional quantum statistics. Q We term g a p as the mutual statistical parameter when a f3 
and self statistical one in the case of a = j3. We call the quantum statistics governed by general values of g a p the 
fractional exclusion statistics. Furthermore, the system is called the generalized ideal gas (GIG) if the total energy is 
given by 



E = Y,N a e c 



(2.4) 



with the one particle energy e Q . Here e a is independent of the distribution {p a } = {N a /G a }. The entropy of GIG 
is given by taking the logarithm of eq. (2.3), 



S = J2 [Pa + N a ) In (D a + N a ) - D a hxD a - N a ln7V Q ] 

a 

= Ga + P<*) 111 + P°^> ~ P a ln P a ~ P*<* hl ^l ' 



(2.5) 



where p* a = D a /G a is the distribution function of dual particle ( hole) of a. p* a is related with p a as 

Pa + 229<xpP0 = i, 



(2.6) 



with 



9aj3 = gapGp/G a , 



(2.7) 



which is obtained by dividing both sides of eq. ( |2.2| ) by G a . 

Thermal average of the particle density in a is given by minimizing the thermodynamic potential f2 = E — TS — 
Y] p, a N a with respect to {p a }. The equation Sfl/dp a = is expressed as 



(l+w a )l[{- 



W/3 



W/3 



9(3* 



exp [(e Va) /T] 



with 



w « = P*a/Po 



(2.8) 



(2.9) 



From eqs. ( p.6[ ) and (2.9), we obtain the thermal distribution p a . In terms of w a , we can rewrite the entropy in a 
simple form as 



S = 22 GaPa [(1 + w a ) In (1 + w a ) — w a lnw a ] 



(2.10) 



The thermodynamic potential has also a simple form: 

n = -T^G a In [1 



(2.11) 



If we consider the case g a p = g a 8 a /3, the above results are further simplified. The particle density is given by 
p a = 1/ (w a + g a ). Here w a is the solution of the following equation: 



exp [(e Q - p a ) /T] = {w a ) 9a (1 + w a ) 



(2.12) 
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III. THERMODYNAMICS OF MULTICOMPONENT SUTHERLAND MODEL 



Thermodynamics of multicomponent Sutherland model was constructed firstly by Sutherland and Shastry. H In 
this section, we summarize the results of thermodynamics of MSM, by following the method in rcf. 9. We consider 
N particle system which consists of Kb component bosons and Kp component fermions. The Hamiltonian is that of 
the modified version of the continuum Sutherland model: frOl 



H = - 



1 y d2 , Z^y A(A - Mjj) 



iKj — - L" \ 6 *3, 



where x t is the coordinate of the i-th particle, L is the linear dimension of the system, and A is a coupling parameter. 
Then My is the coordinate exchange operator which exchan ges the spatial coordinates of i and j-th particles with 
the internal degrees of freedom fixed. For the Hamiltonian (3.1), the energy spectrum in the thermodynamic limit 
has been obtained ||[ll| as 

1 r°° A f 00 f 00 n 2 \ 2 d 3 

E/L=— dkk 2 v{k) + — dk dk'\k-k'\v{k)v{k')+ ________ ^ (3 2 ) 

47Tj-oo 8tt7_ 00 J_ oq 6 

where v{k) is the momentum distribution function and d — N/L = (1/2tt) f^° oo i l (k)dk. Both d and v(k) have 
contributions from each species 

1 f°° 

d = ^d a = — / dk^v a (k), (3.3) 

where v a {k) is the momentum distribution function of particle a, and d a is the a-th component particle density. Since 
each microscopic state is identified with the momentum set occupied by particles, the entropy is given by the same 
form as the free particle case 

S/L = J2^, (3.4) 

a 

with 

±1 f°° 

s a = — dfc[(l±i/ Q )m(l±i/ Q )=F^ a ln^ Q ] • (3.5) 

^ J-oo 

Here and in the following, the upper sign is for boson (a G B) and the lower one for fermion (a € F). Thermal 
equilibrium distribution function is obtained by minimizing the thermodynamic potential 

n = E -TS - L^2fj, a d a , (3.6) 

a 

with respect to each v a . In eq. d3.6| ), we have introduced the chemical potential \x a for each a. The resultant thermal 
distribution functions are 

u a (k) - {cxp [(e(fc) - fi a ) /T] T l}" 1 , (3.7) 

with the one particle energy 

2vr SE k 2 A f°° „.,, ,„ ,,,, TT 2 X 2 d 2 , n _. 

e(k) = T — = - + - dk' \k - k>\ v{k') + —. (3.8) 



Equations ( |3.6| ), (3.7), and ( [3.^ ) give the complete thermodynamics. For later convenience, we introduce the rapidity 
p as 

p(Jfe) = = k + £ f X dfc'sgn (fc - k') v(k'). (3.9) 

ok 2 J_ 00 

Further differentiation of eq. (|3.9|) with k gives 
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d P (k) 
dk 



= l + Xv{k). 



By multiplying both sides of eq. (3.10) by p{k) and integrating them over k, we obtain 

„2 , 



V 



'/2 = 6TT\J2^ [i±^J 



c , 



(3.10) 



(3.11) 



where cq is a constant which shou ld be determined. Since e — > (fc + irXd) 2 /2 and p ^ k + irXd in the limit fc — > oo, we 
obtain cq = 0. By using eq. ( |3 .10 ) and p(k = ±oo) = ±oo, we can rewrite the density for a as 



2tt 



dp 



1 + Xp 



Similarly, the entropy eq. (3.5) is rewritten as 



±1 



dp 



[(1 ± v a ) In (1 ± u a ) T v a In i/ a ] 



2^ 7-oo 1 + ^ 

The internal energy is also rewritten in terms of the integral over p 

1 f°° 
E/L = — dp- 

Derivation of eq. ( |3. 14 ) will be given in Appendix A 
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p A v 
+ Xv ' 



(3.12) 



(3.13) 



(3.14) 



IV. EQUIVALENCE TO THE GENERALIZED IDEAL GAS IN THERMODYNAMICS 

In this section, we show that the thermodynamics of the MSM is equivalent to that of the free particle with FES 
at all temperatures. Firstly we introduce w a as 



I/fa for a £ B 
l/v a - 1 for a e F. 



We rewrite eq. ( |3.11 ) with eq. ( |4.l| ) as 



2 /2 = e + TA^ln 



W a + 1 



By subtracting [i a from both sides of eq. (4.2) and taking the exponential of them, we obtain 



exp [{p 2 /2 - fia) /T) = (1 + w a ) Y[ 



Wp + 1 



with the "statistical parameter " gp a , 



90a 



S a p + X for a £ F 
A otherwise. 



(4.1) 



(4.2) 



(4.3) 



(4.4) 



The quantity p a = + Av) satisfies eq. ( |2 .6|) with eq. ( p.9[) . Comparing eqs. ( |4.3| ) and (|2.8|), we note that p a is 

nothing but the distribution function of free particles with the statistical parameter g a p. 

Next we rewrite the thermodynamic quantities in terms of w a . The total energy is given by 2nE/L = 
E Q S-ao dp (p 2 /2) p a - The entropy of species a is 



2ws a = ± 



dp 



1 + Xv 



[(1 ± v a ) In (1 + v a ) T lni/ a ] 



dpp a [(1 + w a ) In (1 + w a ) — w a \aw a ] 

/oo 
d P [(P* a + Pa) In (p* + Pa) - Pa In p Q - p* lnp*] . 
-oo 



oc- 
oo 



(4.5) 
(4.6) 
(4.7) 
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By using eq. ( |4.7| ), we rewrite the thermodynamic potential as 
Cl/L = E/L — T^~^ s a — [i a d a 



1 f°° 

— I d PPa [p 2 / 2 ~ Va-T(l + w a ) In (1 + w a ) + Tw a lnw a ] 



We then use the relation (4.3) to obtain 

T 

T 

27T l-r 



/CO 
dpp a w a In (l + w^ 1 ) + ^ ffq/3 In ( 1 

/oo 
dp Y, ^ (1 
-oo a 

T f°° / 







(4.8) 



(4.9) 



From the second to the third line in cq. (4.9 ), we have used eqs. ( |2.6| ) and ( pj| ) with G Q = G/3. Compare eqs. (|4.3|). 
(4.7), and ( |4.9] ) with eqs. (2.8), (2.5), and~( 2.11 ). Our results in this section show that the thermodynamics of the 
MSM is equivalent to that of free fractional particle with the statistical parameter g a p and G a — 1 for all a. If we 
set Kb = 1 and Kp — 0, the above results reduce to the thermodynamics of g-on with g — A as expected from results 
in ref. 5. 



V. G-ON DESCRIPTION IN LOW TEMPERATURE REGION 

In low temperature region, we have a simpler picture of the thermodynamics of MSM. In this region, the mutual 
statistics can be absorbed into the self statistics and the system can be regarded as an assembly of g-ons. Before 
getting into details of calculation, first we consider the zero temperature case which will be helpful to understand 
how g-on emerges in low temperature region. For simplicity, consider a simple case of Kb = 1 and Kp = 1 . At zero 
temperature, the distribution functions vb for boson and v-p for fermion are given by 



vn(k) = 27rd B (5(fc), up(k) = 6 (k F - \k\) , 
where fcp = nd-p and 9(k) is the step function: 



9(k) = 



for k > 
otherwise. 



Substitution of eq. (5.1) into eqs. (3.< 

(l + A)fc 2 /2 



and (|3.9|) gives 



e(k) 



7rAd B |fc|+7r 2 A(Ad 2 - 
k 2 /2 + n\d\k\ +TT 2 X 2 d 2 /2 



4) 



/2 



for \k\ < kp 
for &f < \k\ , 



(5.1) 



(5.2) 



(5.3) 



p(k) 



(1 + A)fc + 7rAd B sgn(fc) 
k + nXdsgn(k) 



for I A; I < kp 
for kp < \k\ . 



(5.4) 



From eqs. (p.3|) and (5.4), we obtain the explicit form of e as a function of p 



(X 2 d 2 



Xd 2 ) /2 



7r 2 A (Xd + dp) 2 ^j /2(1 
p 2 /2 



for \p\ <p B 
A) for p B < p < pf 
for pp < \p\ , 



(5.5) 



where pb — TtXdB and pp — tt (dp + Xd). Figure 1 shows the p dependence of e at zero temperature. From eq. ( |3.12 ), 
we can regard p B (F) = 1/ b(f)/(1 + Af) as the rapidity distribution function of boson (fermion). In Fi g. 2 , the rapidity 
dependences of pe and pp at zero temperature are shown by dashed and dotted lines, respectively. II4] The fermion 
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distribution function pp changes at |p| = pp and \p\ = pb while pe changes only at \p\ — pb- We also show the 
profiles of pb{p) and pf{p) at low temperatures quantitatively by dotted-dotted-dashed line and dotted-dashed one, 
respectively. If T is low enough, the temperature dependence of thermodynamic quantities comes mainly from the 
two regions of p: the vicinities of pb and pp. The contribution to the thermodynamics from each rapidity region 
can be expressed in terms of thermodynamics of g-ons. The two kinds of g-ons are independent of each other, if the 
following condition holds: 

Pb(p = Pf)^0 and p F (p = p B ) ~ 1. (5.6) 



The condition eq. (5.6) turns into 

A BF = exp [(/i B - A«f) /T] < 1, (5.7) 
where pb and pp are the chemical potentials of bosons an d fc rmions, respectively. In this case, we define the "low 



temperature region" as the region where the condition eq. (5/7) holds. 

Now come back to the case of general value of Kp. In this case, however, the definition of the "low temperature 
region " depends on the distribution of chemical potentials of each species. In the rest of this section, we consider 
only two cases: the case (1) where all species of fermion have the same chemical potential and the case (2) where 
all species of fermion have different chemical potentials. Generalization to other cases is straightforward. 



A. the case (1) 



In this case, equation (|4.3|) becomes 

exp[(p 2 /2-^B) /T] 



i/'f. 



Up 



(5.8) 



Now we try to reduce eq. (5.8) under the low temperature condition (5.7). First we eliminate wp in eq. (5.8). With 
the help of eq. (4.1), we obtain the following relation: 



w F = (1 + w B ) A 



BF- 



(5.9) 



Substitution of eq. (5.9) into eq. (5.> 



gives 



exp [(p 2 /2 - pb) /T] 



W B^BF 



(1 + wb) X{1 ' Kf} ^ (1 + wbAbf + A BF ) XKf 



(1 + 7« b ) A(1 " A ' fKi (1 + wbA bf ) XKf ' 
Here we define the rapidity region Ho as the region where the following condition is satisfied: 

w b A B f < V- 



(5.10) 



(5.11) 



where rj is a constant which is much smaller than unity but much larger than Abf- In IIo, the right hand side of eq. 
( 5.10 ) turns into 



w b (1 + u;b) 



l-A(l-_fsT F ) * XKp 
^BF • 



From eqs. (fTTol) and (PI) , we obtain 

exp [(e - Co) IT] ~ w 9 B ° (1 + wb) 1 ' 90 for p e n , 
where e = P 2 /2M , M = 1 + XK F , go = A/(l + \K F ), and 

Co = Mb - AK F /i F / (1 + XKp). 



(5.12) 



(5.13) 



(5.14) 



Equation (5.13) corresponds to eq. ( 2.12| ) with energy £oj chemical potential Co ; and statistics go. Here we introduce 
wo as the solution of the following equation: 



G 



exp [(e - Co) /T] = w 9 ° (1 + wo) 1 " 90 for all p. 



(5.15) 



And we r ewrite the boson parts of entropy, density, and energy with Wq. First consider the boson part of the entropy 
(eq. ( |4.5| ) with the upper sign). Since the contribution from IIo (the complementary region of Ilo) is infinitesimal (see 
Appendix B), we obtain 



2ttsb = 



dp 



dp 



[(1 + vb) In (1 + uq) - i/ B In vq] 



(5.16) 



n 1 + ^ 



[(1 + v B ) ln(l + v B ) - v B lnt/ B ] 



dp 



Go 

Gq 



n 1 + XKf + Xvb 
dp 



[(1 + vb) In (1 + i>b) - vb lnf B ] 



n w b + 5o 
dp 

n w o + 5o 



[(1 + wb) In (1 + wb) — wb In wb] 
[(1 + w ) In (1 + w ) - w In w ] , 



(5.17) 



with Go = (1 + X K-p) . It will turn out that Go is the density of states in the absence of other particles. For the 
integrand in eq. (5.17), we can show that 



[(1 + w ) ln(l + w ) - w lnw ] / (w + g a ) ~ for p g n , 



(5.18) 



i n the same way as that in Appendix B. Hence we can extend the integral region to (—00,00) in the last line in eq. 
(5.17) and express 2itsb as 



2nsB — G Q 



dp 



[(1 + w )ln(l + w ) - w In w ] 



00 w o + 9o 

Go / dp [(po + po) In (p + po) - Po hipo - Po m Po] 



(5.19) 



with po = ( w o + go) 1 j which is nothing but the rapidity distribution function of g-on with the statistical parameter 
i?o- Here pj = 1 — goPo is the rapidity distribution function of dual particle (hole) of po. Figure 3 (a) shows the profile 
of po and Pq at zero and low temperatures. Both po and Po change appreciably near the generalized fermi surface 
Po = V 2 -Mo Co- Next we consider the boson density de, which is rewritten as 



2-KdT 



dp 



Go 
Go 



1 + A^ 
dp 



n„ wb + go 
dp 

n w + g ' 



(5.20) 



The integrand (wo + go) 1 in eq. (5.20) is infinitesimal since wo > f?A B p 3> 1 in no. Thus we can extend the integral 

dp 

7o 

(5.21) 



region from Ho to (—00, 00) in the last line of eq. (5.2C) 

2nd B ~ G 
— Gq 



w + go 



dppo 



Similarly, the boson part of the internal energy is given by 

,,2 



2irE- 



B 



dp WTxu)^ G "j_J pp2po/2 - 



(5.22) 
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Next we consider the fermion part. We define the rapidity region 111 as the region where wf > 77 _1 Abf is satisfied. 
Here r\ is a constant which satisfies the inequalities Abf *C f] *C 1. As in the boson case, the main contribution to 
the fermion part of entropy comes only from the region IIi. In III, we obtain 



exp [(ei - fi F ) /T] ~ w^ 1 (1 + rap) 1 91 for p e Hi, 



(5.23) 



with £i = p 2 /2 and g\ = 1 + XK F . Note that this equation is for the g-on with g\ statistics. As in the boson case, we 
introduce w\ as the solution of the following equation: 



exp [(ei - n F ) /T] = wf 1 (1 + wi) 1 91 for all p. 
The fermion part of the entropy has contributions only from IIi and is rewritten as 

dp 



(5.24) 



27TSF = — Kp 
~-Kp 



1 + Xu 
dp 



Ul 1 + XK F u F 
dp 



[(1 — Up) In (1 — up) + u F \nup\ 
[(1 — up) In (1 — up) + up In Up] 



dp 



-oo w i +9i 

OO 



[(1 + w-p) In (1 + Wf) — wf \nwp] 



[(1 + Wi) ln(l + iwi) — tfi Iniui] 



/oo 
dp [Oi + Pi) ln (Pi + Pi) - Pi In Pi - Pi m P*] 
-CXD 



(5.25) 



Here p\ = (wi + 31) 1 and p* = 1 — gipi- The profiles of p\ and p| are shown in Fig. 3(b). The generalized fermi 
surface of p\ is given by p\ — y/2pp. 

In contrast to the case of entropy, the fermion parts of density and energy have contributions from Ho as well as 
from ni. The fermion part of density is 



2irdp = Kt: 



dp 



Up 



1 + Xu 



d P i Kf I d P 

1 + AAp^F 



1 



1 



1 + AAp 1 + XKp + Xu B 

/OO 
dp (pi - XGlp ) . 
-OO 

From eqs. and (Uf), we obtain 

/OO 
dp (GoCoPo + Kp^ppx) . 
-OO 

Similarly, the total energy is given by 

/OO 
dp(G e p a + K F £ipi) ■ 
-OO 

From eqs. ( ^.19 ), (5.25), (5.27), and ( |5.28 ), we obtain the thermodynamic potential as 

/OO 
dp [G In (1 + WQ 1 ) + K F In (l + w^ 1 )] . 
-OO 



(5.26) 



(5.27) 



(5.28) 



(5.29) 



Equation ( 5.29| ) is obtained in the same way as that in the derivation of eq. (4.9). We remark that K F fermi 
components have the same G a = 1 as in the multicomponent FES description, but that the boson component has 
Go < 1 in the g-on description. From eq. ( 5.29| ) , we find that the thermodynamics of MSM is described in terms of 
two kinds of g-ons independent of each other. 
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B. case (2) 



Next consider the case where each species of fermion has a chemical potential different from one another; /ip^ > 
Mf,2 > • • • > Pf,k f > Pb- Here p,F,q is the chemical potential of g-th component of fermion. In this case, we define 
the low temperature region as the region where the following inequalities are satisfied: 



exp [(/i B - Pf : k f ) /T] < 1, 



(5.30) 



exp [(/iF, q +i - MF,q) /T] < 1 for q = 1 ~ K F - 1. 



In this case, equation (O) becomes 

exp [( p 2 /2-hb) /T] = 



(1 + wbY 



q=l 



W F ,q 
1 + W F . 



(5.31) 



(5.32) 



where wp ^ = [fp (e — MF,q)] — 1 and u>b = ^b(e — /^o) From the similar study with the previous case, we find 
that the boson contributions in the present case to the thermodynamic quantities are given by the same expressions 
as those in the previous case, after replacing the expression ( 5.14 ) by Co = Pb — ^Y^, q =i A*F.q/(l + AKp). 



Now we consider the q-th component fermion part. First rewrite the above equation (5.32) in terms of u>f,<j in the 
rapidity region Tl q where wp j9 <C min[A ri9 ] for 1 < r < q < Kf- Here & r ,q = exp [(/iF,r ~ pF.q) /T\. In this rapidity 
region, we have 



1 + WF,r 1 + WF, 9 A 



WF,qA q>r for r < q 
1 for q < r, 



(5.33) 



w B 



1 + i»b 



WF,q\,B - 1 
WF,qAq,B 



Here A q .B = exp [(pF,q — Mb) /T]. Using eqs. ( |5.34 ) and ( 5.33 ), equation ( 5.32 ) becomes 

exp [( e , - C,) /T] ~ (1 + w F ^ 9q for p e n„ 
with e 9 = p 2 /(2M q ), M q = 1 + A (q - 1), Cg = MF, q - AG g X)r=i MF,r, ff g = 1 + AG g , and 

G 9 = AC 1 = (1 + A( (? -1))- 1 . 



(5.34) 



(5.35) 



(5.36) 



Equation (5.35) corresponds to eq. (2.12) with the en ergy e q , chemical potential C, q and statistical parameter g q . At 
the end of this section, we will find that G q in eq. ( 5.36 ) is the density of states of g q -on in the absence of other 
particles. 

We define w q (1 < q < Kf) as the solution of 



exp[(e q -( q )/T]=w^(l + w q ) 1 - 9 « 



for all p. 



(5.37) 



Now we rewrite the entropy, density, and energy in terms of w q . The contribution to the entropy from q-th component 
is given by 



2lIS q = — 



dp 



[(l -,£) In (!-.£)+ if lni£ 



/DO 
dp [(pq + P* q ) In (p q + p* q ) - p q In p q - p* Inp*] 
-co 



(5.38) 



Here p q = (w q + g q ) 1 and p* — 1 — g q p q are the distribution functions of g q -on and its hole, respectively. 
The g-th component fermion density is given by 



2ird q = 



dp 



1 + Xu 



(5.39) 
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Figure 4 shows the profile of the integrand in eq. ( 5.39| ), v q (p) — Vg/(1 + Xv), at zero and low temperature in the 

case where K-p = 2 and q — 1. The generalized fermi surfaces are p q = ^/2M q C sq (q = 0, 1, 2). With the help of Fig. 
4, we see that d q has a positive contribution from the region H q and negative ones from n r (r > q or r = 0) 



2ird q ~ G q 



I — — A V G?/ -* AGq j 

A F 

dpp q - A ^ / dp/5 r - AGq / 

-OO i i «/ — OO J — ( 



dp 



n w + .go 
dpp . 



r= q +l 

Similarly, the energy of q-th component fermion is given by 



2nE q = 



dp 



2 V 

P v q 
2(1 + \u) 



Gn 



AY 



Z r= 9 +l J -' 



dp 



P 2 PB 



From eq. (5.42), we obtain the following result: 



( A F \ A F poo 

?=1 / g=0 J ~ ca 



The total energy is given in a form similar to eq. ( 5.43 ) 

a f 



JV F i> OO 

9=0 " / -°° 



(5.40) 

(5.41) 
(5.42) 

(5.43) 

(5.44) 



From eqs. Q5.37 ), ( 5.38| ), (5.43), and ( |5.44 ), the thermodynamic potential in the low temperature region reduces to 



Af poo 

2%n/L~-T^2G q dplnil + w- 1 ) 

9=0 J' 00 



(5.45) 



VI. DISCUSSION 



In the last section, we obtain the g-on description in the low temperature regime for the two cases ( (1) and (2) 
). Our results imply that the presence of a "magnetic field", which polarizes the system, leads to a g-on description 
different from that obtained in the absence of the magnetic field. As we increase the magnetic field with T fixed, 
thermodynamics at low temperature changes from the one described by the case (1) to another obtained in the case 
(2). Such a statistical crossover of g-on occurs only in the multicomponent model. 

In section 4, we proved that the thermodynamics of MSM is The same as that of particles obeying FES with certain 
statistical parameters. Here we discuss the validity of FES description to other solvable models; recently we proved 
that the thermodynamics of the Haldane-Shastry and long-ranged t-J models are described exactly in terms of the 
thermodynamics of free particle with FES at all temperatures. jl3| Thermodynamics of the Bethe solvable bose gas 
is proved to be equivalent with that of particles obeying the FES in ref. 5. For the isotropic Heisenberg model with 
the nearest neighbor exchange, where the string solutions exist, thermodynamics is described in terms of FES after 
identifying the n string with n-th species (see Appendix C). From the above results, we expect that thermodynamics 
of other Bethe solvable models are interpreted in terms of FES. However, we have not yet obtained the exact FES 
description of the thermodynamics of the multicomponent Haldane-Shastry and long-ranged t — J models. Hence we 
are not sure that thermodynamics of all solvable models are exactly described in terms of FES at the present stage. 

In section 5, we found that g-on description accounts for the low temperature thermodynamics of MSM. On the 
other hand, the validity of g-on description to the lattice versions of the Sutherland model is limited; we found 
that the g-on descriptions provide the correct low temperature properties of Haldane-Shastry, long-ranged t-J, and 
its multicomponent models in the presence of the magnetic field. In the absence of the magnetic field, however, 
thermodynamics of the models are described in terms of the parafermion statistics jl4| , which cannnot be described 
in terms of g-on. Though we have not yet studied the Bethe solvable models at this point, we expect the low 
temperature properties of those models would be described in terms of a picture different from the g-on one. 
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APPENDIX: A 

In this appendix, we derive eq. ( |3.14 ) from eq. ( |3.2[ ). By using eq. ( |3.9| ), we obtain, 

/OO rOO POO 

dkk 2 v + \ / dkku(k) I dk 1 i^(ki)sgn.(k - ki) 
-oo J — oo J —oo 

A 2 r°° 



dkv{k) I dk\v(k{)sga(k — ki) / dfc2^(fc2)sgn(fc — A^) 



The second term in eq. (Al) is rewritten as 

A 



A 



dk / dki (k — ki) sgn (k — ki) v(k)v(k\) 



dk I dki \k — fci| v(k)v(k\). 

— oo J —oo 



The third term in eq. (Al) is rewritten as 

A 2 r°° 



/OO pki 
dki / dk2is(k)is(ki)v(k2)sgn(k — ki)sgn(k — k%) 
-oo J —oo 
2 poo pk rki 

dk f dki / dk2v{k)v(ki)v(k2) 



oo 



2 

A_ 

2 

r 

12 7-oc 

(27rd) 3 A 2 _ 2^ 3 A 2 rf 3 
12 ~~ 3 



oo 



dkv(k) I dk\v(ki) I dk-ivik-i) 



From eqs. flAl]), flA2| ), and ([A3|), we find that 

2ttE/L= I 

J — ( 



p 2 v(k) 



Replacing (ifc — » J*_ dp/(l + Av), we obtain eq. ( 3.14 ) 



(Al) 



(A2) 



(A3) 



(A4) 



APPENDIX: B 



In this App endix, we show that the contribution from the rapidity region IIo to 2ttsb is negligible under the 
condi tion (5.7) in the case (1). In the rapidity region IIo, wb > V^bf holds. The numerator in the integrand in eq. 
(5.16) is rewritten as 



[(1 + Ub) ln(l + ub) - vb hi^B 
w 



= [(l + WB 1 )ln(l + WB 1 )- WB 1 ln WB 1 ] ~0(«;b 1 )<0((A^7 7 )- 1 )«1, (Bl) 
while the denominator 1 + \v is larger than unity. Hence the integrand in eq. ( 5.16| ) is infinitesimal in the region IIo. 
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APPENDIX: C 



Here we consider the one-dimensional isotropic Heisenberg model under the periodic boundary condition: 

JV N 

4 



N N 

H=-J2(<ri-<ri+i-V + Hj2<rt, (CI) 



i=l 



Here er^ is the Pauli matrix at the site i and N is the number of site. Magnetic field H is applied parallel to z-axis. 
The thermodynamics of the model eq. (CI) were first constructed by Takahashi by using the string hypothesis. 



In ref. 15, the thermodynamics is given as follows. The internal energy E is 

E/N = J2 dke%{k)v n {k) - H, (C2) 

i J — oc 



n=l 

with 



where k is the rapidity and v n (k) is the rapidity distribution function of the string with the length n. The entropy S 
is given by 

°° />oo 

s/n = Y, dk [("" + O ln ly* + <) - v » ln - < ln <1 • ( C4 ) 

n=1 J-oo 

Here i/* is the rapidity distribution function of hole of n-string. Among v* and v m , there is the following relation: 

°° />oo 

G n (k)=v*(k)+Y, / dk' 9nm (k-k')u m (k'), (C5) 

m=1 J-oc 

where G n (k) — n/-K(n 2 + k 2 ) and g nm {k) is given by 

ffnmO) = /|n- TO |( fc ) + 2 /|„- m | +2 (fc) + 2/|„_ m | +4 (fc) 4 h 2/ n+TO _ 2 (fc) 4- fn+ni(k), (C6) 

where /„(fc) = n/n(n 2 + k 2 ), for n > and fo(k) — 5(k). Here we introduce p n — v n (k)/G n (k) and p n (k) = 



u*(k)/G n (k). In terms of yO n and p*, the relation (C5) becomes 

1 = -°n+E / dk'g nm (k-k') Pm (k'), (C7) 

m=1 J-oc 

with <jr„ m = g nm G m /G n . Similarly, the energy and entropy are given as 

00 />oo 

E/N = J2 dkG n s° nPn (C8) 

n=1 J -co 

°° />0O 

S / N= H / dfcG « [(-°« + -°«) ln 0« +Pn)- Pn ^ Pn - p* n hi p*] . (C9) 

1 J — oc 



From the expressions (C7), (C8), and (CE), we can regard the thermodynamics of isotropic Heisenberg model as that 



of GIG with energy s n , statistical parameter g nm (k), and the density of states G n . 
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FIG. 1. Rapidity dependence of one particle energy e(p). Two dots represent the points where the derivative with p is 
discontinuous. 



FIG. 2. Rapidity distribution functions of boson ps(p) and fermion Pf(p) are shown in the case where Kb = Kf — 1. The 
dashed and the dotted lines represent pb(p) and Pf(p) at zero temperature, respectively. The dotted-dotted-dashed line and 
dotted-dashed one represent Pb{p) and Pf{p) respectively at low temperature. 



FIG. 3. (a) Rapidity dependences of po{p) and po(p) at low temperature are shown by dashed and dotted lines, respectively, 
(b) Rapidity dependences of pi(p) and pl(p) at low temperature are shown. In both figures, results at zero temperature are 
also shown, for comparison. Note that pi(p) has no structure at p — po in contrast to Pf{p) in Fig. 2. 



FIG. 4. The profile of the integrand in eq. ( 5.3S ). vi(p) in the case of Kf = 2 is shown at low temperature. The generalized 
fermi surfaces for po, pi, and p2 are represented by po, pi, and p2, respectively. For a reference, the result at zero temperature 
is also shown by dotted line. 
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